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Abstract 

We show that the equation of quasigeostrophic (QG) potential vor- 
ticity conservation in geophysical fluid dynamics follows from Hamil- 
ton's principle for stationary variations of an action for geodesic motion 
in the /- plane case or its prolongation in the p - plane case. This im- 
plies a new momentum equation and an associated Kelvin circulation 
theorem for QG motion. We treat the barotropic and two-layer baro- 
clinic cases, as well as the continuously stratified case. 



To appear in Physics of Fluids 

1 Introduction 

The quasigeostrophic (QG) approximation is a basic tool for the analysis 
of meso- and large-scale motion in geophysical and astrophysical fluid dy- 
namics |Q| . Physically, QG theory applies when the motion is nearly in 
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geostrophic balance, i.e., when pressure gradients nearly balance the Cori- 
olis force. Mathematically, in the simplest case of a barotropic fluid in a 
domain T> on the plane M. 2 with coordinates (x±,X2), QG dynamics in the 
/3-plane approximation is expressed by the following evolution equation for 
the stream-function ip of the geostrophic fluid velocity u = z x ~Vip 

where d/dt is the partial time derivative, A is the planar Laplacian, T 
denotes rotational Froude number, la, 61 = Jp2£l i s the Jacobi bracket 
(Jacobian) for functions a and b on M? and is the gradient of the Coriolis 
parameter, /, taken as / = /o + (3x2 in the /3-plane approximation, with 
constants (3 and /q. Neglecting j3 gives the /-plane approximation. The QG 
equation Ql.ip may be derived from the basic equations of rotating shallow 
water flow by a proper nondimensionalization and subsequent asymptotic 
expansion with the Froude number given by the square of the ratio of the 
characteristic scale of the motion to the deformation radius (cf. Pedlosky 
[i] and, for recent discussions, Allen and Holm and Stegner and Zeitlin 
[3]). Alternatively, equation ( |l.l[) may be written in terms of the potential 
vorticity as 

^ + u-Vg = 0, 9 = A^-^ + /. (1.2) 

In this form of QG dynamics, the basic property of potential vorticity con- 
servation on geostrophic fluid parcels becomes obvious. 

Much is known about the mathematical structure of the QG equations 
fll.ip and ( |1.2| ), In particular, they are Hamiltonian with a Lie-Poisson 
bracket given in Weinstein [|]] 

{F,G} = - J d Xl dx 2 q[^, d -^-] , (1.3) 

where jx = q — f '. In terms of the variable jx, the Hamiltonian for QG is 
expressed as 



H = - I dx x dx 2 (\Vip\ 2 + Fi> 2 



^ J dxidx 2 n {T - A)~V + \Ylf i™-dx, (1.4) 



where 7« is the z-th connected component of the boundary. In what follows we 
shall limit ourselves by the cases where T> is either a torus (periodic boundary 
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conditions) or the whole plane M 2 with decaying boundary conditions and, 
thus, the boundary terms may be ignored and 5H/5fi = ift. Consequently, 
the Lie-Poisson bracket ( |1.3| ) gives, after integration by parts, the dynamical 
equation for fi, 

^ = {^H} = -[^ q ] = -u.Vq, (1.5) 
in agreement with the QG potential vorticity equation (|1.2|). Casimirs of the 



Lie-Poisson bracket (L3) are given by C$ = f dx\dx2&(q) for an arbitrary 
function $ and satisfy {C$,H} = for all Hamiltonians H(fi). Level sur- 
faces of the Casimirs C$ define coadjoint orbits of the group of symplectic 
diffeomorphisms of the domain of the flow [Q] , || . 

The QG approximation may be generalized to the baroclinic case of sev- 
eral superimposed fluid layers [0. For the simplest two-layer model, the 
equations of motion describe the potential vorticity conservation in each 
layer 

+ U W • Vg (i) = 0, u« = z x , i = 1, 2 , (1.6) 
where potential vorticities are defined as 

'W = dW(A^W + (-l) < J^ i )(^ 1 )-^ 2 )) + /) * = 1,2, (1.7) 



(f 



and dW and 

are the nondimensional layer depths and rotational Froude 
numbers, respectively. In the two-layer case, these quantities are related by 
= J r ( 2 )rf( 2 ). The symplectic structure of the iV-layer QG equations 
is well known II. The Poisson bracket for multilayer QG is also Lie-Poisson 
of the type ( |1.3[) . The Hamiltonian in this case is again quadratic and equal 
to the sum of the kinetic energies of the layers and the so-called available 
potential energy. In the case of two layers the Hamiltonian is given by 

H = \f d Xl dx 2 ( |VV W | 2 + ^(^ (1) (i {1) +^ (2) d (2) ) h/>«-^ (2) f 

J ^i=l,2 

(1-8) 

There is also a QG approximation for a continuously stratified baroclinic 
fluid In terms of the (positive) "static stability" parameter S(z), the 
QG equation for a continuously stratified fluid is 

d f . , d 1 dtb\ , d 1 dib „\ . 

« ( A * - ssw + u • V ( A * - ifzwd + ; ) = (L9> 
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Again this is potential vorticity conservation on geostrophic fluid parcels. In 
fact, this is the standard equation for oceanic synoptic motions In this 
equation, u = z x V?/> is the geostrophic velocity for a stratified fluid and the 
stream function ip(xi, X2, z) depends on both horizontal and vertical coordi- 
nates in the domain Px [zq, z\}. The Hamiltonian formulation of baroclinic 
continuously stratified QG dynamics with a Lie-Poisson bracket of the 
type is given in Holm Q. The Hamiltonian in the case of appropriate 
boundary conditions in the vertical direction is P] 



1 



H = - dx\dx2 I dz 
2 Jv Jz 



|vv>r + 



2 , 1 W 



S(z) \dz 



(1.10) 



where the term 2 s(z) (§z ) 2 * s * ne density of available potential energy in 
stratified QG theory 

In this paper, we provide new variational principles for QG dynamics 
in several cases, and formulate the corresponding momentum equations and 
Kelvin circulation theorems in these cases. These variational principles for 
QG dynamics are potentially useful in observational data assimilation and 
studies of the transport and mixing properties of QG flows. Our derivations 
are based on the previously mentioned fundamental relation between the 
QG equations and the group of symplectic diffeomorphisms (area-preserving 
mappings) of the domain of flow T>. Related ideas have recently been dis- 
cussed from the Euler-Poincare viewpoint in Holm, Marsden and Ratiu H|. 
Our approach is similar to that of Virasoro || who suggested using the gen- 
erating function for symplectic diffeomorphisms as a generalized canonical 
coordinate with the subsequent construction of a Lagrangian via the Legen- 
dre transform of the known Hamiltonian expressed in terms of this variable. 
However, such a procedure can be practically realized only in case where 
the incompressibility constraint may be explicitly solved, i.e., for diffeomor- 
phisms that are sufficiently close to identity and, thus, describe only the 
initial stage of dynamical evolution. We avoid solving the incompressibil- 
ity constraint and, instead, use the parameters of the group of symplectic 
diffeomorphisms of the domain of the flow as dynamical variables. 

A frequently used variational formulation of two-dimensional incompress- 
ible flows is the Clebsch formulation; see, e.g., Seliger and Whitham [10|, 



Holm and Kupershmidt [O], Marsden and Weinstein 12], Salmon |13j and 
references therein. This formulation can indeed be extended successfully 
to derive the quasigeostrophic equations. However, we decline using the 
Clebsch representation here because it is not universal. That is, not every 
coadjoint orbit of a given group can be represented by a given set of Clebsch 
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variables. For a discussion of this point see, e.g., Zeitlin [14], [15|, |lq] ; a 
more detailed analysis of Clebsch representations for the case of symplectic 
diffeomorphisms will be given elsewhere. 

2 Preliminaries: Variational principles for Euler 
dynamics 



We begin with the case of /-plane dynamics, for which equations (|1.2| ) and 
(1.6) for potential vorticity convection may be interpreted as geodesic equa- 



tions on the group of symplectic diffeomorphisms in the same way as for 



the classical 2d Euler equations discussed in Arnold [17|. The difference 
is in the choice of metric which is determined by the relation between the 
stream- function and the (potential) vorticity |5| , |l8| . Without giving rigor- 
ous differential-geometric and functional-analytic meanings to this interpre- 
tation we shall simply discuss it heuristically, by using the guideline that in 
formulating a geodesic variational principle on a manifold it is convenient to 
deal with two sets of variables: those belonging to the linear (co-)tangent 
space (the dual of the Lie algebra of planar divergenceless vector fields in our 
case); and those belonging to the configuration manifold itself (the Lie group 
of area-preserving mappings in our case). In fluid dynamics, this means us- 
ing both Eulerian and Lagrangian coordinates. See, e.g., Holm, Marsden 
and Ratiu Q for a discussion of the equivalence of these two descriptions 
from the viewpoint of variational principles defined on Lie groups and Lie 
algebras. 



Lagrangian description. As shown in Arnold [17|, the square of the 
geodesic distance on the group of symplectic diffeomorphisms is given by the 
kinetic energy of the fluid, and by explicitly introducing the incompressibility 
constraint with the help of a Lagrange multiplier (the pressure) one obtains 
the action principle for 2d incompressible fluid dynamics 

5S = 0, with S = j dt j dhdl 2 Q|x| 2 +p(J- 1)^ . (2.1) 

In this expression, x = (xi, x^j, the variables x\{l\, I2; t),X2(h, h',t) are con- 
sidered as coordinates of fluid particles marked by Lagrangian labels h,l2, 
and J is the Jacobian of the map (h,l 2 ) —> (x\,x 2 ), with J\ = dx l /dl A , J = 
det(J^). Finally, overdot ( ) denotes time derivative following the fluid 
parcels, i.e., at constant Lagrangian labels. The standard Lagrangian equa- 
tions of motion for an ideal fluid of a constant density are (see, e.g., Lamb 
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1) 

x + -^ = 0, J=l. (2.2) 
ax 

These equations are found from Hamilton's principle upon varying S in x, 
imposing the constraint J = ^fe^y = 1 (area preservation) on the map 
(h,h) (%i,%2) and using the known properties of Jacobian determinants 

d(a,b) d(a,b) d(x 1 ,x 2 ) , d(a,x 2 ) da 

and a /_ ITT = "5Z~ - ( 2 - 3 ) 



d(h,l 2 ) d(xi,x 2 ) d(h,l 2 ) d(xi,x 2 ) dxi 



The action 5 in equation ( |2.lD is quadratic, since the pressure term vanishes 
on the constraint manifold. This is appropriate for a geodesic principle. 
(For a discussion of this point, cf. f20|| .) As is well-known, equations ( [2,2j ) 
transform to the standard Euler equations for the velocity field u = x upon 
using the formula 



<it \dt 

Note that the Lagrangian conservation of vorticity (Kelvin circulation theo- 
rem in 2d) follows from the general structure of the equations ( |2.2D . Elimi- 
nating the pressure in ( |2.2D gives 

^1 _ = o (2 5) 

dxi dx 2 ' 

which is equivalent, in view of (|2.3|), to u = 0, where to = §f^ — §fj = ^-curlu 
is the vertical component of vorticity. 

Eulerian description. An alternative way to treat the action principle 
( |2.lD is to consider it from the very beginning in the Eulerian (x\,x 2 ) plane, 
as 

S = jdtL, L = J dxidx 2 ^D\u\ 2 - p(D - 1) , (2.6) 

and to perform variations at fixed xi,x 2 and t of the inverse maps (x±, x 2 ) — *■ 
(h,l 2 ), viewed as Eulerian fields l A (x\,x 2 ,t), A = 1,2, which determine u 
and D by (summing on repeated indices) 

xiA aiA 

% F = -u i Dt Dt = l- D = det(Df). (2.7) 
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Variations with respect to l A and p of any action whose Lagrangian is of the 
form L{u,D;p) then yield |IJ, ||, 



5S 



dt I dxdy { D{D~ 1 ) i A 5l A 



d 1 8L 1 5L j /5L\ 
dtDfru 1 + D5ul U ^ ~ \6D~),i 



5p(D - 1) 



(2. 



where d/dt = d/dt + u • V is the material derivative. Stationarity of the 
action S(u,D;p) under variations in l A also results naturally in Kelvin's 



circulation theorem for any action of this form, since |21], [22] 

•(fx = 



d r i 8L , 

— <b — — -dx. 



dt ./ 7 ( t ) D 5\i 



7(<) 



d 1 SL 1 SL „ , 

1 :VW J 

IdtDSu D8w> 



i{t) 



0, 
(2.9) 

where 7(f) is any closed loop moving with the fluid velocity u. For the action 
( |2.6[) we have the variational formulas 



1, 



— u 

2 1 1 



P- 



(2.10) 



1 SL 5L 

Dlhl~ u, Id 

Hence, from the vanishing of the coefficient of 5l A in ( [2.8| ) we find the Euler 
equations for the velocity field, 







— u — u x curlu + VI p -\ — |u 



1 



0, 



(2.11) 



(2.12) 



m ' " V ' 2' 

upon using the fundamental vector identity of fluid dynamics, 

(b • V)a + aj Vb> = (curia) x b + V(a • b) , 

with, in this case, a = and b = u. 

Of course, geometrical considerations are not essential in formulating this 
variational principle, because x(l, t) represents the fluid particle positions 
and equation ( |2.2D is just Newton's second law with a constraint. In fact, 
this variational principle was first written by Lagrange himself [23]. However, 
geometrical considerations are sometimes helpful in its interpretation. For 
example, 2d incompressible fluid dynamics in the presence of the spatially 
dependent Coriolis force (fo+(3x2)uxz may be obtained from this variational 
principle by a "prolongation" |@] ; ||, pi] ] of the momentum density (which 
is equal to the velocity in the present context) in either (^l|), or (^), i.e., 
by a shift 

5L 5L , . 

► h D~R, or, equivalently in this case, u — > u + R, (2.13) 

ou ou 
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where 

K = y(xi + Px\x 2 ), so z ■ curlR = f{x\,x 2 ) = fo + @X2- ( 2 -14) 

Geometrically, this shift corresponds to a central extension of the Lie algebra 
of symplectic diffeomorphisms 0. The prolonged action may be rendered 
quadratic by introducing an auxiliary gauge field using the Kaluza-Klein 
construction Q. 



3 Variational principles for QG dynamics 

As in the case of classical Euler dynamics described in the previous section, 
we shall use the parameters of the symplectic diffeomorphisms of the flow 
domain as basic dynamical variables for QG. Thus, at this level of approxima- 
tion we introduce "quasigeostrophic particles" which move at the geostrophic 
velocity u = z x Vt/> and follow the quasigeostrophic fluid trajectories, Fin 
pre which result from a certain filtering of the primitive equations. These 
QG fluid trajectories will be described as functions of Lagrangian coordi- 
nates given by xi(li, l 2 ; t), x 2 (h, l 2 ', t) in V. Such quasigeostrophic particles 
are commonly used in the context of another geophysical fluid dynamics 



approximation, namely semigeostrophy [24] 



Lagrangian representation. To apply the geometric approach of the 
previous section to QG dynamics, we consider a generalized quadratic action 
on the /-plane of the form (cf. equation (2.1) 

dt I dhdi 2 di' x di' 2 \Ut{hM) ■ G{*{hM)Mi'iA)WiA) 



+ dt dhdl 2 p(J - 1} 



(3.1) 



where £7(x, x') is the kernel (Green's function) of a bounded self-adjoint 
scalar integral operator. Varying this action with respect to x at fixed 1 and 
t gives the following (Eulerian) equations of motion 



d_ 
dt 



dx'idx 2 Q(x, x')x' 



dx' 1 dx' 2 -^ : G(-K,x')x' i )xi + 



X 



dp 



0, (3.2) 



where the change of variables from 1 to x is understood. As before in (|2.5| ), by 
eliminating the pressure and using identities fl2.3j ) we obtain a conservation 
law on parcels, q g = 0, where 



1g 



d 
dx 2 



dx' 1 dx' 2 Q(x,x')x'i + 



d 
dx\ 



dx' 1 dx' 2 Q(x, x')x' 2 , (3.3) 
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This expression recovers the potential vorticity for QG in the f-plane, i.e., 
qf = Aip — Tip in equation (|1.2| ), when Q is chosen to represent 1 — .FA -1 
on the domain T> and x = z x V?/> is used. The passage to the /3-plane case 
may be accomplished by the prolongation ( |2.14 ). 



Eulerian representation. Prolonginging 5L/5u — > 5L/5u + DR, in the 
action (|3.1| ) written in the Eulerian velocity representation with integral 
operator Q = (1 — JFA -1 ) yields S = J dt L, with Lagrangian 



L 



dx\dx2 



-Dm • (1 - TA' l )u + Du ■ K(x 1 , x 2 ) - p(D - 1) 



(3.4) 



and produces the following variations at fixed x±,x 2 , and t, 



1 5L 

D~5u~ 
SL_ 

Id 

SL 
5p 



R + u 

1 

2 

-D + l. 



A _1 u 



A- L {Du) 



2 2D 
u- (1 -^A _1 )u + u-R-j?. 



(3.5) 



Hence, from equation (2.8) for action principles of this type and the funda- 
mental vector identity ( 2.12| ), we find the Eulerian QG motion equation, 



— (l-^A- 1 )u-uxcurl((l-^A- 1 )u+R)+V^+-u-(l-^A- 1 )uJ =0. 

(3.6) 

where we have used the symmetry of the Laplacian and substituted the 
constraint D = 1, imposed by varying p. The curl of this equation yields 



dt 



+ u-V( At + /) + (/i + /)V-u = ! 



(3.7) 



where fj, = z ■ curl(l — ^-"A _1 )u. The constraint D = 1 implies V • u = 
(from the kinematic relation dD /dt + 'V- Du = 0) and when u = z x is 
substituted, the equation for n = Aip — Tip yields the QG potential vorticity 
equation Ql.2| ). The QG version of the Kelvin circulation theorem follows 
from formula (12.91) as 



d_ 

lit 



I ( (1 - FA'^u + R) ■ dx = , 

/ 7 (t) v ' 



(3.8) 



for any closed loop ^{t) moving with the QG fluid velocity u under the 
Eulerian motion equation ( |3.6| ). The QG motion equation ( 3.6| ) is derived in 
Euler-Poincare form via the Kaluza-Klein construction in 



D. D. Holm and V. Zeitlin — Hamilton's principle for QG 



10 



For a geodesic stationary principle, the Lagrangian and the Hamiltonian 
coincide |2C]. In our case, Legendre transforming the Lagrangian (3.4) gives 
the Hamiltonian 



H 



dx\dxo —Du 
12 



1 - JFA-^u + p(D - 1) 



(3.9) 



Thus, the Lagrangian and the Hamiltonian do coincide for QG, in the f- 
plane theory and on the constraint manifold D = 1. Taking u = z x V-0 
and using the divergence theorem yields the Hamiltonian for QG in terms 
of the stream function, 



H 



D=l 



— I dx\dx2 
v 



(3.10) 



thus reproducing the QG Hamiltonian flL4J) for the Lie-Poisson bracket for- 
mulation discussed in the introduction. 



The multi-layer QG equations. A similar variational program may be 
carried out in the case of the multi-layer QG equations. In what follows we 
shall confine our attention to the two-layer case, the further extension to N 
layers being straightforward. In order to build an action principle, we should 
take into account that, geometrically, two-layer QG dynamics corresponds 
to geodesic flow on the direct product of two groups of symplectic diffeomor- 
phisms uM , each group being related to the motion in one layer. We, thus 
introduce two sets of Lagrangian fluid trajectories x and y, the pressures 
p( x >y) i n each layer and the corresponding integral operators 

g(xx) = _ ^Wa" 1 ), g(w) = d {y) (l - JF^A" 1 ), 

g{xy) = hjr{x) d {x) + jriy) ^ A -\ ( 3-11 ) 



Here the superscripts x,y correspond to 1,2 in (y^) , (|L7|) . The two-layer 
analog of the action (|3.l| ) then may be written as 



s = f dt f d/ 1 d/ 2 <^[ix(/ 1 ,/ 2 )-^)(x(/ 1 ,/ 2 ),x(r 1 ,^))x(/;,^) + 

^(hM)-Q (xy \AhM)Mi\A)m'iA) + 
\y{hM-Q {yy \y(hM)Mi[A)m[A) 

+ j dt j dhdl 2 (p {x) {J {x) - l)+p iv) (J iy) - 1)). (3.12) 
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By performing the variations and eliminating the pressure as before we find 
potential vorticity conservation in each layer. For example, in the x layer 
we have 

^W = 0, (3.13) 
where (using = 1 and = 1) 

q(x) = + u = / ^4s ( ** ) (x,x'K + / M^^WM- 

(3.14) 

A similar expression holds with obvious changes in the y layer. It is easy to 
check that both these expressions coincide with the definitions Ql.7| ). Thus, 
two-layer QG dynamics is obtained from Hamilton's principle for the action 
( |3.12| ). Again, one could start by transforming the action ( 3.12j ) into the 



Eulerian velocity representation by using two sets of Lagrangian labels and 
arrive at the same multi-layer QG equations. 

Continous vertical stratification. Consider an action of the following 
Eulerian form in three dimensions (x\, X2,z) with vertical coordinate z in 
the interval [zo,zi], 



S = I dt dx\dx2 \ dz 

IV Jz 



r ^Du-(l-C(z)A- 1 )u+Du-K(x 1 ,x 2 )-p(D-l) 



(3.15) 

Let the operator C(z) depend only on z, so C(z) commutes with the planar 
Laplacian, A. Hamilton's principle for this action yields a motion equation 
corresponding to (|3.6|) 



— (1 - £(z)A- x )u - u x curl((l - £(z)A- x )u + R 

+ v(p+-u-(l-£(^)A- 1 )u) =0, (3.16) 



with V = (d/ 8x1,8/8x2), as before. Thus, by equation Q2.9D , the Kelvin 
theorem for stratified QG dynamics has the same form as the single layer 
Kelvin theorem fl3,8| ) on each horizontal surface parameterized by the level 
z = const. Namely, cf. equation (|3.8|), 

— I f(l-£(z)A~ 1 )u + RVdx = 0, (3.17) 
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for each closed loop j(t) moving with the horizontal QG fluid velocity u, gov- 
erned by the Eulerian motion equation ( |3.16| ). This corresponds to "pseudo- 
potential vorticity" conservation. The curl of the stratified Eulerian QG 
motion equation fl3.16|) yields, cf. equation ( |3.7|) , 



d_ 

dt 



Aip-£(z)ip) + u- V [Aip~ £{z)ip + f) = 0, 



(3.18) 



where u = z x Vip(xi,X2, z) is consistent with D = 1 and the operator £(z) 
is still to be specified. Choosing 



£(z) 



d 1 d 



(3.19) 



dz S{z) dz' 

which contains the (positive) "static stability" parameter S(z), produces 

~d 2 ip 



d dip d dip d 
dt dx\ 8x2 8x2 8x\ 



d 2 ip | #V | d / 1 dip\ + ' 
8x\ 8x\ dz \S{z) dz J 



0. (3.20) 



This is the standard equation for oceanic synoptic motions, cf . Pedlosky Q , 
p.364. 

Legendre transforming the Lagrangian in (3.15) yields the (constrained) 
Hamiltonian, 



H 



dx\dx2dz 
dx\dx2dz 
dx\dx2dz 



u ~ - -Du • (1 - £(z)A -1 )u — Du ■ R + p(D - 1) 
du 2 



^Du ■ (l - £(z)A-^u + p(D - 1) 
^D\Vip\ 2 - ^DVip • £(z)A- x VVj + p(D - i; 



(3.21) 



When this Hamiltonian is evaluated on the constraint manifold D = 1, it 
becomes the conserved energy, 



E = — I r/./'j d.vodz 



\Vip\ 2 - ip£(z)ip 



(3.22) 



after integrating by parts in (x%, X2) and discarding the horizontal boundary 
term. For the choice of the operator £ in ( 3,19j ), we then find 



E 



— I dx\dx2dz 



— I dx\dx2 



itt/i2 1 f 9 ^ 

ip dip 1 " 



2 1 



S(z) dz 



(3.23) 
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which is the conserved energy for ideal synoptic scale motions and the Hamil- 
tonian for the Lie-Poisson formulation of baroclinic QG theory, cf. equation 
(p. .ICQ and Holm J?]]. A choice of C(z) extending that in equation ( |3.19| ) 
yields the modified QG theory of White |25[|, namely 



£(2) = (.I + B )^)(|- B )-^ (3 - 24> 



where the definitions of the symbols may be obtained in [25] . For a discussion 



of the Lie-Poisson Hamiltonian formulation of White's modified QG theory, 



see Holm [26|. 



4 Discussion 

There are several potential applications of the QG variational principles 
presented here. For example, by studying the second variation of the action 
and, hence, the separation of the geodesies (cf. [pC[] ) one may attack the 
problem of Lagrangian stability of QG flows and study dispersion of the 
quasigeostrophic particles which may be considered as a test of accuracy and 
predictability of QG flows. Also, in the context of numerical simulations, 
one may try to find proper discretizations directly at the level of the action, 
as proposed in Holm et al. (See Marsden and Wendlandt f28fl , [29] for 



recent advances in this direction.) Finally, as QG dynamics is frequently 
used in assimilating atmospheric and oceanic data, the variational principles 
presented here could be useful in this domain, e.g., if optimization procedures 
were applied directly to the QG action, cf. Bennett |pCf| . To facilitate their 
utility, we have expressed these QG variational principles in both Eulerian 
and Lagrangian coordinates. 
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